Production of resonances in a thermal model: 
invariant- mass spectra and balance functions 
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Abstract. We present a calculation of the 7r + 7r - invariant- mass correlations and the 
VO \ pion balance functions in the single-freeze-out model. A satisfactory agreement with 

the data for Au+Au collisions is found. 
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Thermal models of particle production in ultra-relativistic heavy-ion collisions 
turned out to be very successful in describing hadron yields [1-4], spectra [5,6], the 
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elliptic-flow coefficients v 2 [7,8], and the HBT radii [7,9]. In view of this fact it is 
challenging to study correlation observables in the thermal approach and compare the 
results to the data. Particular examples of such observables are the 7r + -7r - invariant- 
mass spectra [10] and the balance functions [15,16]. The latter describe correlation 
between opposite-charge particles in the rapidity space and are closely related to charge 
! fluctuations [17-19]. 

In order to analyze the 7r + 7r~ invariant mass spectra measured in Ref. [10] with the 
like-sign-subtraction technique, we make the assumption that all the correlated pion 
pairs are obtained from the decays of neutral resonances: p, Kg, rj, r/, fo/cr, and / 2 . 
The phase-shift formula for the volume density of a n-n resonance with spin degeneracy 
g is given by the formula [13] 
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which was used in Ref. [11], as well as by Pratt and Bauer [14] in a more recent 
analysis. We have performed a calculation in the framework of the single-freeze-out 
model of Ref. [5], including the flow, experimental kinematic cuts, and decays of higher 
resonances. The full results are shown in Ref. [11]. In the left panel of Fig. 1 we 
compare the model predictions to the STAR data. The mass of the p meson was scaled 
down by 10% in accordance to the experimental evidence for the dropping mass [10,20]. 
Our results were filtered with the detector efficiency correction (we are grateful here to 
P. Fachini). We can conclude that the model does a very good job in reproducing the 
gross features of the experimental invariant-mass spectra. In the right panel of Fig. 1 we 
show the model predictions for the tranverse momentum spectra for several resonances, 
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Figure 1. Left: single-freeze-out model vs. the data of Ref. [10] for the 7r + 7r~ 
invariant- mass correlations. The model calculation includes the decays of higher 
resonances, the flow, the kinematic cuts, and the detector efficiency. Right: parameter- 
free model predictions and data for the transverse-momentum spectra of several 
resonances. 

confronted to data. The overall agreement is impressive. Note that all the parameters of 
the model (two thermal and two geometric) had been fixed with the help of the spectra 
of pions, kaons, and protons, such that there was no more freedom left in the analysis 
of resonances. 

The balance functions, measured by the STAR Collaboration [21,22], are defined 
by the formula 



B(5,Y) = - 



1 \ (N + -{S)) - (N ++ (6)) , (N. + (S)) - (N._(S)) 
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where N+_(5) denotes the number of the opposite-charge pairs such that both members 
of the pair fall into the rapidity window Y with relative rapidity \y 2 — yi\ = 5, N + is the 
number of positive particles in the interval Y , etc. The measurement [21] showed that 
the widths of the balance functions are smaller than expected from models discussed 
in Ref. [15]. This problem was discussed by Bialas in Ref. [23] in the framework of the 
quark coalescence model [24]. Here we present the results of the calculation of the 7r + 7r~ 
balance function [25] based on the model or Ref. [5]. In a more recent paper Pratt et 
al. [26] showed that the measured balance function may also be satisfactorily reproduced 
in the blast-wave model. 

In our approach the tt + tt~ balance function has two contributions related to two 
different mechanisms of the creation of an opposite-charge pair. The first one (resonance 
contribution) comes from decays of neutral hadronic resonances, whereas the second one 
(non-resonance contribution) is related to other possible correlations. We assume that 
the second mechanism forces the two opposite-charge pions to be produced at the same 
space-time point and with thermal velocities. The first contribution includes neutral 
resonances which have a 7r + 7r~ pair in the final state (we explicitly include Ks, rj, rj 1 \ p°, 
a;, a, and / ); here the correlations are completely determined by the kinematics. The 
details of the model and other technical remarks are given in Ref. [25]. 

According to the above discussion, the 7r + 7r~ balance function can be constructed 
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Figure 2. Balance functions for the pions in the thermal model calculated for four 
different centrality classes and compared to the experimental data of Ref. [21]. The 
normalization of the model curves was adjusted in each case and is listed near the plot 
labels. 

as a sum of the two terms, 

B{6) = B K {6) + B mi {6). (3) 

The functions B R (5) and £> NR (5) resulting from our model calculation are presented in 
Ref. [25]. The value of the temperature used in the calculation was T = 165 MeV, and 
the expansion parameters of the model had been fitted to the spectra of hadrons [5] . This 
procedure yields the average transverse flow of 0.5c. One can observe that the widths 
of the two contributions are similar. The calculated total width, (5) = 0.66, turns out 
to be somewhat larger than the experimental value for the most central collisions. The 
STAR result for the most central events (c = - 10%) is (S) = 0.594 ± 0.019, for the 
mid-central (c = 10 - 40%) (S) = 0.622 ± 0.020, for the mid-peripheric (c = 40 - 70%) 
(5) = 0.633 ± 0.024, and for the peripheric (c = 70 - 96%) (5) = 0.664 ± .029. Such 
dependence of the width of the balance function on centrality cannot be reproduced 
in our model via changing the transverse flow within the limits consistent with the 
single-particle spectra. 

In Fig. 2 our results are compared to the experiment. The normalization of the 
model curves was adjusted in each case, since we were not able to take into account, in 
a more sophisticated way, the effect of the limited detector efficiency and acceptance. 
On the other hand, the kinematic cuts in pseudorapidity and transverse momentum 
were included [25] . The shapes of the model balance functions agree well with the data 
except for the most central case where the theoretical width is slightly larger. The dips 
of the experimental balance functions at very small values of S are caused by the HBT 
correlations, not included in our approach. We also note that the effects of the detector 
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efficiency may influence the width of the balance function [26]. 

In conclusion, we state that the single-freeze-out model gives a satisfactory 
description of the pion invariant-mass correlations and of the balance functions. The 
results for the transverse-momentum spectra of resonances are in remarkable overall 
agreement with the experiment. 

We are grateful to Patricia Fachini for numerous helpful discussions and for 
preparing the figures. This research has been supported in part by the Polish State 
Committee for Scientific Research grant 2 P03B 059 25. 
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